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C D ' Abstract. 4-dini intrinsic (material) Riemannian metric G of the material 

C D 4-D space-time continuum P is utilized as the characteristic of the aging pro- 

^Nj , cesses developing in the material. Manifested through variation of basic mate- 

t, • rial characteristics such as density, moduli of elasticity, yield stress, strength, 

O ' and toughness., the aging process is modeled as the evolution of the metric G 

^^ ^ (most importantly of its time component Gqo) of the material space-time P 

^^ , embedded into 4-D Newtonian space-time with Euclidean metric. 

VP i The evolutional equation for metric G is derived by the classical variational 

approach. Construction of a Lagrangian for an aging elastic media and the 
derivation of a system of coupled elastostatic and aging equations constitute 
the central part of the work. The external and internal balance laws associated 

^ ■ with symmetries of material and physical space-time geometries are briefly re- 

^+ ' viewed from a new viewpoint presented in the paper. Examples of the stress 

f^ ' relaxation and creep of a homogeneous rod, cold drawing, and chemical degra- 

f**^ , dation in a tubing are discussed. 

o 
\o 
o 

1. Introduction 

Qh' We seek to develop a model of inelastic processes in the aging materials by 

(-H , employing a 4-D inner material metric tensor G as the aging (damage) parameter 

of a material continuum. Aging here implies any variation in the chemical make-up, 
i.e., chemical degradation, phase transformation, phase coarsening, nucleation, and 
growth of micro-defects such as dislocations and voids, shear bands, crazes, micro- 
cracks, etc. Material engineering and failure analysis indicate that, in addition 
to the stress and strain tensors, a parameter of state (the "aging" parameter) 
^ , is needed to represent on a continuum level the sub-micro and micro-structural 

j^ ■ changes of material. A kinetic equation for the evolution of the aging parameter 

will represent the aging process of a material. The equations of evolution for the 
material metric G are the Euler-Lagrange equations resulting from a Variational 
Principle. The conjugate force of the evolution of metric G (and of the related 
quantities characterizing the properties of the material) is the Energy-Momentum 
Tensor of Elasticity introduced by J. Eshelby (see Sec. 7 below). 

A 3-D material metric g has long been employed as an internal variable in con- 
tinuum mechanics. For example, it was used for studying the duality of material 
and physical relations of the Doyle-Erikson type in article ^ , the thermodynamics 
of a continuum in [2j , and in where the curvature of material metric g defined 
by a uniformity mapping of a uniform material was employed as the driving force 
of the material evolution. We use the 4-D material metric G as an additional state 
parameter that reflects the aging process. G is introduce with the largest covariance 
group allowed by the condition that a small vicinity of each point of the material 
preserves its topology during the aging process (see Sec. 2 below). We consider the 
3-D material metric g on the slices Bt of constant physical time as one of the main 
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dynamical variables following the ADM (Arnovitt, Deser, Misner) presentation of 
G (see ^E] or Sec. 2. 3 below). In that respect, we follow the tradition of the cited 
works. What is new in our work is that the smaller (in comparison to the General 
Relativity) covariance group of the Lagrangian allows us to use the lapse function 
S and the shift vector N as independent dynamical variables reflecting the proper 
material time scale and the intrinsic material flows respectively. 

This aging parameter is justified by the observation of shrinkage associated with 
aging and the subsequent material density variation as well as a change of the res- 
onance atomic frequencies and characteristic relaxation times measured in macro- 
scopical studies. In other words, the internal length and time scales change with 
aging when compared to the corresponding absolute (physical) scales. The most 
sensitive indicator of aging is a variation of an intrinsic material time scale. The 
measurement of time in the laboratory as well as in material (intrinsic time) can 
be accomplished by several methods, the most common of which is the use of os- 
cillating processes such as those found in clocks with a pendulum or crystal-based 
timepieces. Another way of measuring time is the use of a unidirectional evolution 
of state. In medieval Europe, for example, time was measured by burning a candle 
which had numbered and colored beeswax strips. Still another method is associated 
with relaxation processes which require an excitation input to enable a fading re- 
sponse. Electronic relaxation generators employing the discharge of a capacitor and 
the fading luminescence of phosphorus are both examples of relaxation processes, 
which are well suited for measuring intrinsic time scale changes because they reflect 
atomic or interatomic events. 

Consider an external excitation of a material which responds with a specific 
change in its state. The decay or fading of the response constitutes the relaxation 
process. The decay can be described by an exponential function (within certain 
limits) e'' "^^ where t is time and tq is the time constant characterizing the rate of 
relaxation. Usually tq becomes smaller with an increase in temperature or decrease 
in pressure. Phosphorous fades more slowly at colder temperatures, for example. 

In section 2 we discuss the kinematics of a media with a variable Riemannian 
metric G in a 4-D material space-time P, embedded into 4-D Absolute (Newton's) 
space-time Af^ with the Euclidean metric H . Thus, the 3-D "ground state" metric 
tensor is introduced together with the proper time lapse function and the mate- 
rial shift vector field. We consider mass conservation law in section 3. Elastic and 
inelastic strain tensors -E^' ; E™ are introduced in section 4 as a measure of deforma- 
tion and the "unstrained state" respectively. The Lagrangian describing inelastic 
and elastic processes in media is discussed in section 5. A variational formulation of 
aging theory and the Euler-Lagrange equations (equations of elasticity coupled with 
the aging equations) are considered in section 6. We present the combined system 
of elasticity and aging equations in section 7 and discuss special cases of the aging 
equations in section 8. Corresponding to the material and laboratory symmetries, 
we consider the space and material balance laws in section 9. In section 10 the 
Energy-Momentum balance Law and the decomposition of the Energy- Momentum 
tensor into components, including the Eshelby tensor and terms related to the ag- 
ing processes are presented. In the final section we explore the application of this 
model to the basic inelastic processes- unconstrained aging, stress relaxation, and 
creep of a homogeneous rod. 
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2. 4D KINEMATICS OF MEDIA WITH A VARIABLE METRIC. 

In this section we introduce the basic elements of the kinematics of a continuum 
with a variable metric, including material space-time P, 4D material metric G, 
4D deformations t/i, slicing of the material space-time by the surfaces of constant 
physical time i?0.f, and total, elastic and irreversible strain tensors. 

2.1. Physical and Material Space-Time. Let us consider the 4-D Euclidean 
vector space M = R x R^ (physical space-time) with the standard Euclidean 
metric H . There exists the volume form d^v corresponding to this metric. 

We select global coordinates x%i = 1,2,3, in the physical space R'^ and x^ = t 
on the time axes R. We have H — dt^ + h = dt^ + ^^ dx"^ ^ . 

Hyperplanes t = c are endowed with the 3D Euclidean metric h induced by H. 
We extend 3D tensor h to the degenerate (0,2)-tensor h in M, taking fiQi — hiQ = 
hoQ = 0. 

A solid is considered here, in a conventional way as a 3D manifold B with the 
boundary dB, i.e. a set of "idealized" material points. Wc will use local coordinates 
X^ , I = 1,2,3 which, incidently, may be global coordinates induced by a reference 
configuration i.e., a diffeomorphic embedding (j)o 3 : B 1 — > R"^ (,6 ). 

Cyhnder P = R x i? (with the coordinates (X° = T, X^, / = 1, 2, 3)) is equipped 
with the 4D Riemannian metric G (material metric) with the components Gij 
relative to the coordinates X° — T, X' , I = 1,2, 3. Space {P, G) is further referred 
to as the material "space-time". 

Metric G defines the 4D volume form dV — \J\G\d'^X, where \G\ is the determi- 
nant of the matrix (Gij). 

An example of such a material metric G can be constructed as follows. Ex- 
tend the reference configuration c/jq 3 to the diffeomorphic embedding 0o : P ' — *■ 
M, (j>o{T,X') = (T,a;* = c/)}^ ^iX')). Let Go be the metric (t)oiH) (here and below 
we denote by (j)*Q the puUback of a covariant tensor Q by the differentiable mapping 
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(p). In the coordinates (X,T) the matrix of the metric Gq is , . Denote 

V L) nij J 

by dVo = y^\Go\d'^X = (/)3((iu^) the 4D-volume element defined by the metric Gq. 

Projection tt : P 1 — > B along T-axes plays the same role in the construction 

below as in the relativistic elasticity theory (E1,|H])- In particular, we require in- 

variance of Lagrangian theory with respect to the automorphisms of the bundle 

(P, TT, B) (diffeomorphisms of material space-time P onto itself, projecting to B, so 

that material points retains their identity during the material evolution) preserving 

the direction of the flow of the "intrinsic" time (see below), but not with regard to 

the whole group of diffeomorphisms of P as in Gravity Theory. 

2.2. Deformation History. The history of the deformation of the body B 
is represented by a diffeomorphic embedding (j) ; P 1 — > M of the material 
space-time P into the physical (Newtonian) space-time M (see Fig.l). 

A deformation history </> for which t = 4>'^{X) = T will be called "synchro- 
nized" . The synchronization can practically performed for relatively slow defor- 
mation processes (in comparison to sound wave velocity). 

Using the deformation cj), we introduce the slicing of the material space-time P 
by the level surfaces of the zeroth component of (j) 

B^ = {B^^t = cfP -\t)} = {{T,X) e P\4>\T,X) - t}. (2.1) 
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Figure 1 . 4D Deformation History and the Field of Flow Vector 



For a synchronized deformation Bcf,^t — {{T,X) e P\T = t}. 

There is a time flow vector field W0 in P, associated with the slicing B^,t of 
the space-time P (gj, [3)- This (future directed) vector field represent the flow of 
"intrinsic" (proper in Relativity Theory) time in the material. Lifting the index in 
the 1-form d^" with the help of the metric G and normalizing obtained vector field, 
we define the time flow vector field as 

"* = II ,,nii ■ (2-2) 



The norm of the f-form d(jf is defined as Hd^^lP — {G^^cjfj^cjfgY/'^ (summa- 
tion agreement by repeating indices is used here). Thus, u^ is the unit vector 
G-orthogonal to the slices i?0,t. In the local coordinates X^ , 



G^^fs 



d 



yG^V>^ 



dX^' 

B 



For the synchronized deformations does not depend on 

G^o d 



'-G 



(2.3) 



(2.4) 



Additionally, if the metric G has the block-diagonal form in the coordinates (X° = 
T,X^) (shortly, BD - metric), we have Uq = [Goo]^^-^- 

Let Mo = ^ be the flow vector associated with the metric Go and the corre- 
sponding 3D shcing Bq. 

We require fulfillment of the following condition ensuring the irreversibility of 
the flow of time: 

<u^,uo>G >0- (2.5) 

Deformation history cf) for which the condition (2.5) is satisfied is called admissible. 
In coordinates (X^) this condition reduces to the following simple inequality 

^% > (2.6) 

and, therefore is a restriction on the deformation history only. 

For a synchronized deformation history (j), this condition is trivially satisfied. 
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Time component (jP of the deformation history may be excluded from the hst 
of dynamical variables by an appropriate " gauging" . Namely, we use the invari- 
ance of Lagrangian under the automorphisms of the bundle (P, tt, B) to make the 
deformation history synchronized. 

An automorphism F : P i — > P, X^ = F^ (Y^) of the bundle (P, tt, B) determines 
the diffeomorphism of the base B that can be considered as a change of variables 
XI = F^{Y^). 

In the new variables, the condition (2.6) takes the form qyo = 4'^iPfi = ^xi-^.o > 
(since F^ , I — 1,2,3 do not depend on F°). Thus, the class of admissible 
deformation histories (f> is stable under the action of the subgroup Auf^ (P) of all 
automorphisms of P with Pg > 0. 

The group Aut'^{P) of automorphisms of the bundle contains two subgroups. 
One is the subgroup TC of the "time change" proper gauge diffeomorphisms (X" — 
T;Xi,I = 1,2,3) I — > {F{T,;X'',J = l,2,3),X'f,/ = 1,2,3) for an arbitrary 
smooth function F{Xi) with Fq > 0- The other one (denoted D[B)) consists of 
the lifts to the slices P^.t of the manifold P of the orientation preserving diffeomor- 
phisms of the base B (group of such transformations of P is denoted Dif f^{B)). To 
lift a diffeomorphism we use (diffeomorphic) projections tt^^j = 7r|B_^ ^ : P^ f i — > B. 
If (p is synchronized, lifted diffeomorphisms do not depend on T . 

Any automorphism of the bundle (j) e Aut^ (P) generates the time independent 
diffeomorphism 0b of the base B, that is element oi Diff^{B). Lifting this element 
to the element of D{B) we represent the group Aut~^{P) as the semi-direct product 
of the normal subgroup TC{P) and the subgroup D{B). Thus, we have proved the 
first of two following statements 

(1) Automorphisms group Aut^{P) is the semidirect product of the subgroup 
D{B) ^ Diff^{B) of orientation preserving diffeomorphisms of the base 
B and the subgroup TC of proper gauge transformations ^_f of the fibers 

Cf : {{X° = T; X^ 7 = 1,2,3)^ (FiT, ; X\ J = 1, 2, 3), X' , 1 = 1,2, 3)) 

with Po > 0. 

(2) For any admissible history of deformations one can choose a transforma- 
tion ^p £ TC such that the history of deformation (j) ° £,f is synchronized. 

To prove the second statement let be an admissible history of deformation. De- 
fine the element P G TC{P) as follows: P : (T, X^, / = 1, 2, 3) i — > {(t>°{T, X^ , J = 
1, 2, 3), X^, I = 1,2,3). Then, = 01 o P where (pi is another admissible history of 
deformation with the same components 0*,i = 1,2,3 and the identity component 
(t>1{T,Xi,I = 1,2,3) = T. Transformation P is admissible since P"o = 0°o > 0, 
thus P g TC{P). Apparently, the deformation history 0i is synchronized. 

Therefore we may restrict our consideration to the synchronized histories of 
deformation keeping in mind that the covariance group of the theory reduces from 
the group Aut~^{P) to the group D{B) of time-independent orientation preserving 
diffeomorphisms of base B. 

2.3. ADM-decomposition of Material Metric, Lapse and Shift. Slicing 
B^t, generates the (3,l)-decomposition of the material metric G employed (for 
a Lorentz type metric) in General Relativity (@],|S]). Specifically, the sandwich 
structure of the part of the manifold P bounded by the surfaces B^^t and B^^t+At 
allows one to introduce the time-dependent lapse function TV and the shift vector 



ALEXANDER CHUDNOVSKY AND SERGE PRESTON 



To+AT 




Cp,Tn+AT) 



^B 



(<p,T:ny 



Figure 2. Lapse Function S and Shift Vector N. 



field N tangent to the shces B^.f such that the metric is block-diagonahzed in the 
moving coframe {dT,dX^ + N^dT): 



ds^ = GijdX'dX' = gu{dX' + N'dT){dX'^ + N'^dT) + S'^dT^ 



(2.7) 



Matrix representation of the material metric tensor G and inverse tensor G ^ 
have in these notations, the forms 



Goo 
G/o 



Go,/ 
Gij 



NaN^ + S^ 

Nr 



Nj 

gij 



G™ 
G^o 






1 
_n' 



TV' 
-,IJ I A'-^A'-' 



g- + 



S2 



and g 



(2.8) 
is the correspond- 



where g is the 3D-metric induced by G on the slices B^ 
ing inverse tensor. In these notations vJGI = S'-\/|5|. 

In what follows we assume that the 4D-deformation history 
Thus slices B^j has the form T = t — const. 

In these notations, the flow vector ug and the corresponding 1-form have the 
form 



is synchronized. 



UG^^dT-—dxA, u^a^SdT. (2.9) 

The last formula gives the " material time differential" dr = SdT for the material 
metric G. In our context, the coordinate X dependence of lapse function S(T,X) 
accounts for heterogeneity of material aging in different points of the solid. On the 
Fig. 3 the local observer at different points of body at different moments of time T 
sees the different rate of the local time in comparison with the laboratory clocks. 

Moreover, the lapse function S can be considered as an intrinsic measure of 
material age, associated with its cohesiveness. It can be normalized to be equal 
1 in the reference state of the solid. As a result of energy dissipation in various 
inelastic processes, material loses its cohesiveness with aging. In the formalism 
presented here it is manifested in slowing down of the material (intrinsic) time, i.e. 
increasing of S{X, T). 

Here we do not consider thermodynamics. However, monotonic increase of 
S{X,T) resembles and can be linked to the principle of non-negative entropy pro- 
duction of the thermodynamics of irreversible processes: ^||^ > 0. 
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Figure 3. Illustration of Time Rate S (Lapse Function) by Vari- 
ation of Time Interval (r — tq) = S{Ti — Tq). 




Figure 4. Propagation of Phase Boundary Represented by Shift 
Vector N 



The requirement S > 1 leads to the strong constraints on the form of the " ground 
state term" of the material Lagrangian L^, (see Sec. 5). 

In this context, the shift vector field N in the metric G reflects a propagation of 
the phase transition or chemical transformation boundary through the material as 
reflected, for instance in the mass conservation law (see below). 

The separation of the evolution of 3D material metric g, material transformation 
process in Bt, characterized by N and the inhomogeneity of the local time dr = 
SdT, are the main reason for introduction (3+1) ADT-representation of 4D-metric 
in Gravity Theory (0]). In addition, an adoption of this view leads to a very clear 
separation of the "physical" degrees of freedom in the canonical formalism and to 
the explicit hyperbolic formulations of Einstein Equations ( JOIE])- 



3. Mass conservation law 

The mass form dAI = podV deflned in P is introduced here, in addition to 
the volume form dVc of metric G. The reference mass density po, defined by this 
representation, satisfies the mass conservation law ([11,13) 



£u,dM = d{iu,dM) ^ 0. 



(3.1) 



Here Cu^ is the Lie (substantial) derivative of the exterior 4-form dM in the direc- 
tion of the vector fleld u^. Recall that the Lie derivative of a differential form ui 
along a vector filed u is deflned as £„a; = ^(/•(WJi^o, where (p^uj is the puUback of 



the form uj by the flow t 



— exp{tu) of the vector field m ([HI)- 
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Equation (3.1) is equivalent to the condition divcipau^) = 0, where divergence 
is taken with respect to the volume form dV . In local coordinates the Mass Con- 
servation Law has the form 




= 0. (3.2) 

Due to the properties of the metric G and the deformation </>, the material space- 
time manifold P is foliated by the phase curves of the flow vector field u^ and thus 
the value of the reference mass density yOo(0, X) at T = uniquely defines its values 
for all T > 0. 

In the synchronized case 0° = T, ||d(/)0|| = VG^ = S'~\ G*^° = "^^ G"" = 
S~^ and (3.2) take the form of the following balance law 



(Pov^).o = E O^^POV^ , ■ (3.3) 



From (3.3) we note that the shift vector field N can describe the matter (den- 
sity) flow due to the some internal processes such as the phase or chemical 
transformations. 

If, in addition to being synchronized, the material metric G is also in the BD- 
form {N = 0), the flow term in (3.3) disappears and the mass conservation law is 
equivalent to the following representation of the reference mass density in terms of 
its initial value po(0, X) — po'- 



Mr,A-',^M0,.v,^^.-^, (3., 



where G(0, X) is assumed to be Euclidean metric. If metric G does not changes 

dpo 
dT 



with time T, we get the classical local mass conservation law -^ = 0. 



4. Elastic, Inelastic and Total Strain Tensors 

In this section we introduce the principal quantities characterizing both elastic 
and inelastic deformation processes. Total deformation is presented as a composi- 
tion of elastic and inelastic ones and is integrable. Its elastic and inelastic "com- 
ponents" are non-integrable, in general, but might be such in special situations 
(see Sec. 11). We recall that the presentation of total deformation as a composition 
of this type was studied in different forms in many works (O^l^J, to name a 
few). What is new here is the 4D-approach to this decomposition and direct defi- 
nition of elastic, inelastic and total strain tensors in terms of material metric gt as 
an independent dynamical variable, reference (undeformed) Cauchy metric go ^-nd 
the current Cauchy metric C^{(1>) rather then using the "deformation gradients" 
(integrable or not) of elastic and inelastic (plastic) deformations. 

NaN"^ Nj 



Slicing Bifi _t of P defines the covariant tensor 7 = G~Ufi,®u^ — , 

(|Z|j@])- Here and later the sign s over = means that this equality is true in syn- 
chronized case. Tensor 7 induces the time dependent 3D-metric gt on the slices 
-B^ ,t (see, for example, [TT]V 
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To obtain the expression for gt in material coordinates X^ , we notice that the 
tangent vectors 

6---Tfaxo+ax., /-1,2,3 (4.1) 

form the basis of the tangent spaces to the shces B^i. In this basis, gt is given by 

QAB = 9{£,a,^b) ^ Gab -GnB-T) ^^o — tt + ^oo ' , n ? ' A,B = l,2,i. (4.2) 

<^'o '?^'o '/'o 

For a synchronized deformation 7 = G — {G^^)^^dT ® dT and gt is simply the 
restriction of 4D-metric G to the slices Bt=c, i-e. gt ij = G/j|T=t, see (2.8). 

Denote by ht the 3D metric on the leaves B^t induced by the metric Gq (that 
is by the tensor jg = Gq — uq iS) uq) . 

Associated with the tensor 7 there is the projector ((l,l)-tensor) 

r<10 / r, Q \ 

n = G-'^ = I-u;®u^=^ I-—dxi®dT^i^^j "/^j. (4.3) 

on the tangent spaces to the slices B^^t, last equality being true for synchronized 
deformations 0. 

Let us consider the puUback 6*4(0) = (f)*h oi the degenerate tensor h by the 4D- 
deformation mapping (f). Tensor C4(0) is degenerate in P, its kernel is generated 
by the vector (j)~^(-S^). In coordinates (X-^) we have 






C4W1J =1 ^ ^, ,j ^3 J, :, ^j I • (4.4) 



The spatial part of this tensor is the conventional Cauchy-Green strain tensor C3((/>) 
of the Elasticity Theory. Components of this tensor with indices (OJ) and (10), 
/, J = 1,2,3 have the form velocity x deformation covector (see JHl)- (00)- 
component of ^4(0) is the square of the material velocity V = (03^(^) =^ qT'Sj^- 

4.1. Elastic Strain Tensor. Here we define the 4D (l,l)-elastic strain tensor 
in P. Wc will do it first in linear approximation and then, using logarithm of a 
(l,l)-tensor function, in another way, more suitable for large deformations. 
We start with the following, conventional definition: 



1 

2 


/ 5-2||V||2-^-27V/(V,0-,), 


5-2(V,0-,)„-5-27V^(0-^0-^), \ 


-S-^N'\\Y\\l-N'+ 


-5-2iV^(V,0,,)^+/^C3(0)A'j- . 


V+(5-2iV^7V^+g^^)(0-^,V)„ 


-,55-5-2iV^iV^(0;^0-,)^ ) 






(4.5) 



This tensor contains the square of material velocity and the shift vector field. 
Having in mind the general, dynamical situation it is more appropriate to use the 
following tensor as the proper Elastic Strain Tensor 

^"' (0) ^H^-^'H^ iHG-VG4(0)-7))n='' i ( IK I J. ?\ i,Tl IKr< ,1 xl 

(4.6) 
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Here 11 is the projector on the sUces i?0,f defined in (4.3). The last equahty is vahd 
in the synchronized case. Notice that the basic invariants Tr{A'^) for the tensor 
(4.6) are the same as for the tensor Cj,{(f)) — g. 

For the simphcity we use the same symbol E"^^ for the restriction of this tensor 
to the slices B^t. 

Tensor E{(f)Y^ is a measure of the deviation of the Cauchy metric C^{(j)) of the 
actual state from the "ground state" G. For a synchronized deformation (f) and a 
material metric G with the zero shift vector, E"^^ = i(C(0) — gt) has the form of 
the conventional elastic strain tensor. 

Remark 1. The deformation (p is essentially 3-dimensional in the sense that 
only the spatial Euclidean metric qq — h 'u\ B is deformed. The 4D-tensor Gi{(f)) = 
(ffh defines the degenerate metric in the material space-time P. It is instructive 
to compare G/i{(j)) with the (degenerate) tensor 7 = G ~ u^ ^ u^. The elastic 
strain tensor E'^^ measures the deviation of C4{(f>) from 7 on the slices -B^,*- Thus, 
the scheme presented here is essentially different from relativistic elasticity theory 
([Z|;[H]) as well as from 4D version of conventional elasticity theory. 

We see from (4.6) that E'^'^ = if and only if the following two conditions are 
fulfilled: 

fl) 9ij = ^3(0)/,/ = (t)lih)ij, 

\2)iV^=5^^(0:o,0:^). ^ ^^ 

In particular, metric g coincide with the Cauchy metric induced by deformation 
(p and is flat. 

If E'^'- = then E'^'- = if and only if in addition to the conditions (4.9) the 
following conditio is fulfilled 

3) 11^11'=. 9'^(0;o>;7>(0;o,'^;,/)- (4-8) 

4.2. Inelastic Strain Tensor. Now we introduce the inelastic strain tensor in 
linear approximation 

(last equality being true in synchronized case) and the total strain tensor E*°* of 
the body at each given moment T to characterize the deviation of the deformed 
Euclidean metric (/)*/i|b^j from the initial (Euclidean) 3D-metric h [h being the 
restriction of Go to the slices i?0,t) 

^*°* = iG-i(C4(0)-7o). (4.10) 

Tensor ij*°* can be represented as the sum of the elastic strain tensor E"^^ and of 
inelastic strain tensor i?'": 

&°^ = E"^ + E"^ . (4.11) 

To obtain the corresponding decomposition for the 3D strain tensors we apply 
projector 11 to the total and inelastic strain tensors. In particular we introduce 

E" ^ UE'-U -■\{^,_ /.,^^^s ,, _ ;»•,,,,,, ) , (4^12) 
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As a result, we get from (4.11) the corresponding decomposition of "3D total strain 
tensor" 

n£;*°*n = e"^ + s™. (4.13) 

Restriction of these tensors on the 3D slices i?0_t leads to the more conventional 
(i-dependent) version of this decomposition. 

For a synchronized deformation history, restriction of E^" to the slices Bt takes 
the form 

E-'Ib, = lg-\g - go) - ^(/ - g~'go), (4.14) 

that describes the decline of 3D material metric g from its initial (reference) value 

go = 00 'i- 

Another way to define strain tensors, more suitable for description of large de- 
formation is to take 

E'^' = hniG-'Ciim^E^'' = lln{G^'G)),&"' = hniG^'C^m). (4.15) 

We can define E'^'', E^", i?*°* correspondingly, using projector H. Strain tensors, 
defined in such a way will, in some simple cases, enjoy the same additive relations 
as (4. 11), (4. 13). On the other case, if elastic deformation happens in the directions 
different from the principal axes of inelastic deformation, relation between these 
deformations becomes more complex. 

The relationship between these definitions and those of the linear approximation 
above is established by using the fact that for a couple A, B of (0,2)-tensors such 
that A is invertible, ln{A^^B) « A^"^ [A— B) provided A—B is small enough. Thus, 
when linear approximation is allowable, first definition is the good approximation 
of the second. For instance 

ln{g-'Csm^ln{I + g-\C,{q^)-g))^g-\C3{<l>)-g), (4.16) 

provided C3 (0) — 5 is small. 

4.3. Strain Rate Tensor. One can also define the material elastic strain rate 
tensor as follows 

E = Cu,E'\ (4.17) 

as well as inelastic strain rate tensor 

E =/:„,i?". (4.18) 

In the case where G = Gq and 0° = T, elastic strain rate tensor defined in (4.17) 
has, the same spatial components as the conventional strain rate tensor ([Sj)- 

Denote by G the the Lie derivative G = Cu^G of the metric tensor G with 
respect to the flow vector u^. Then the calculation of the Lie derivative in (4.17) 
results in the following relation 

e' = Cu,E-' = -G-^GE-' + iG-i(C4'(0) - K), (4.19) 

where K — Cu^l is the extrinsic curvature tensor of the slices -B^,*. 

Remark 2. Here we are using material coordinates and tensors only. In order to 
obtain the corresponding "laboratory" quantities (seen by an external observer), 
one defines the laboratory (Euler) Elastic Strain Tensor 

e^-' * = 0^^^t1 ^Ei ^ (4.20) 
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Figure 5. Decomposition of Deformation History. 

and recalculate all the other quantities accordingly. 

Figure 5 presents the above decomposition of total deformation into the inelastic 
and elastic deformations. 

The actual state under the load at any given moment T results from both elastic 
(with the variable elastic moduli) and inelastic (irreversible) deformations. The 
"ground state" of the body is characterized by the 3D-metric gt- This state is the 
background to which the elastic deformation is added to reach the actual state 

(in)- 

Transition from the reference state to the "ground state" that manifests in the 
evolution of the (initial) Euclidean metric h to the metric gt cannot be described, 
in general, by any point transformation. Transition from the "ground state" to the 
actual state at the moment t also is not compatible in general. Yet the transition 
from the reference state to the actual state is represented by a diffeomorphism (^t. 

Here we are considering the material 4D-metric G and the deformation 
(or elastic strain tensor E*^' {(/))) to be the dynamical variables of the field theory. 
The reference mass density po is found (for the synchronized deformation (j) and the 
BD-metric G) by the formula (3.4) if its initial value po{T = 0) is known. In this 
study we consider mainly the quasi-static version of the theory, i.e. inertia forces 
and kinetic energy are assumed to be negligible. 

5. Parameters of Material Evolution, metric Lagrangian 

In examining the processes of deformation and aging of a solid with a synchro- 
nized deformation history (j) we use both general (Gij) and ADM {S, N, g) notations 
for the 4D material metric G. 

Following the framework of Classical Field Theory (^3) we take a Lagrangian 
density C{G, (p) referred to the volume form dV as a function of 4D-matcrial 
metric G, its invariants (with respect to the group of Dif f'^{B) of the orientation 
preserving diffeomorphisms of the base B, see above) and the elastic strain tensor: 

L(G,£;^')- 

The Lagrangian L(G, E''^)dV is represented as a sum of the two parts: the metric 
part Lm{G) and, as a perturbation of the ground state, the elastic part Le{G,E^'') 
associated with elastic deformation 

L = L^{G)+L,{E''',G). (5.1) 
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Metric Lagrangian Lm{G) in (5.1) is introduced to account for the "cohesive 
energy" or strength of the solid state, the strain energy of "residual strain" and the 
energy of the change associated with a evolution of material properties in time, for 
instance material aging processes of phase transitions. 

Metric Lagrangian Lm is the sum of several terms with the coefficients that may 
depend on the 3D volume factor \gt\ (actually |5t|/|5o|) and the lapse function S. 
These volume factors are associated with the solid state ability to retain its intrinsic 
topology in contrast to the fluid and gaseous states. 

First term of Lm is the "ground state energy" F{E™,S, ||iV||^) (shortly GS) - 
initial ("cohesive") energy (per unit volume). 

The second (kinetic) term in L„i (see (5.4) below) is the function of invariants 
of the tensor K = CugI of extrinsic curvature of the slices Bt in the material 
space-time P ( piTT]'!. 

In the ADM notations, the (l,l)-tensor K has the following form 

where (*) represents the terms which do not enter the invariants of K and £, ■ g = 
-Mpg — Cj^g. Lie derivative C^g of the tensor g with respect to the vector field N 
is calculated on each 3D slice Bt for a fixed T. 

In the case of a block-diagonal metric G (no shift: N — 0), ^ — i 

therefore, K is, essentially, the time derivative of the metric gt'. 

Therefore, tensor K represents the rate of change of intrinsic length scales that 
reflects the aging processes. It shall be noted that K is also related to the elastic 
strain rate (see (4.19)). 



Lm{G) also may depend on the shift vector field N through its norm ||A^||g — 
NaN^ (entering the "ground state energy" F) and, possibly, divergence diVg{N) 
and the "proper time derivative" CuN. 

We also include the term reflecting the residual strain energy (incompatibility) 
of -E*°* which is accounted for by the scalar curvature R{g) of 3D material metric 

9- 

Summarizing the above assumptions we construct the metric Lagrangian Lm as 
the scalar function of the parameters listed above: 

Lm = FiS, \\N\\% E^ + x{K) + adiVg{Nf -f /3i?(gO- (5.4) 

Function x of invariants of tensor K^j (" dissipative potential", comp.^21) corre- 
sponds to the energy of inelastic processes in the material °. 
Coefficients a,/9 may depend on S and \g\. 



"Initially f ll7l ') we've considered function x to be a quadratic function of invariants 
Tr{K), Tr{K^) but as the examples of stress relaxation and creep in a rod demonstrate 
this function should be chosen differently, corresponding to the material studied. In partic- 
ular, if the Dorn relation f] = D{exjr'^ — 1) between the stress a and the strain rate f] 
{ri{t) is the volume preserving part of inelastic strain ) is to be obtained, one should take 
XW = cx+ fjjln{%) - i(l + %)ln{l + %). 



14 ALEXANDER CHUDNOVSKY AND SERGE PRESTON 

In the case of a homogeneous media or m ID case, the scalar curvature R{gt) of 
the metric gt is zero and the corresponding term in (5.4) vanishes. 

Given the diversity of the material properties and the different conditions (of 
loading, boundary, forces, heat, etc.) of inelastic processes affecting the material 
it is especially important to choose the material Lagrangian of different materials 
appropriately. It appears as if the different conditions activate different "layers" 
of structural changes for a given material and, correspondingly, turn on terms in 
the "ground energy" and in the "dissipative potential" that are responsible for 
the given type of aging. For example, the slow process of unconstrained aging in 
a homogeneous rod (see sec.ll. or J^l) is overcome by the scale processes of a 
stress relaxation or creep each of which begins in a different loading situation after 
the strain energy (density) reaches a (different) activation level. For these two 
processes both ground energy F and the dissipative potential x(^) have the same 
form different from those for unconstrained aging. 

Remark 3. The scalar curvature R{G) of the 4D metric G can be expressed, by 
the Gauss equation, as the combination of scalar curvature of 3D metric g and of 
invariants of its extrinsic curvature: R{G) = —{tr{K'^) — (trK)'^) + R{gt), up to 
a divergence term ( ,11, ) . As a result, the above form of Lagrangian for an aging 
media (5.4) is a generalization of the Hilbert-Einstein Lagrangian R{G)^y\G\ of the 
General Relativity (0]). By breaking of the invariance group of general relativity to 
the smaller group of automorphisms of the bundle P —^ B we can use more general 
form of metric Lagrangian. 

The perturbation of Lagrangian due to elastic deformation is taken in the form 
of the Lagrangian of Classical elasticity (j^. Sec. 5. 4) 

L^iE-'.G) = ^\\V\\l - p^f{E^\ G) - p^U o <^, (5.5) 

where Poll^llh = PoJ2ij ^jj'^!o'^^o is the density of kinetic energy, / is the strain 
energy per unit of mass, U is the potential of the body forces. Strain energy / is 
assumed to be a function of two first invariants of the (l,l)-strain tensor _E^'. Strain 
energy may depend on the metric G through the invariants of g^ g^, S, vector field 
N, scalar curvature R{g) etc. 

Because we are considering a quasistatic synchronized theory here we ignore the 
inertia effects and, therefore, omit the kinetic energy term in (5.5). 

The strain energy density in linear elasticity is conventionally presented as follows 

where /U, A are the initial values of the Lame constants f|18|). 

We assume that Strain Energy / and the "ground state" term F are independent 
of each other. Yet, in Appendix A we introduce a scheme where elastic deformation 
(elastic strain tensor i?*^') is considered as (small) perturbation of (large) inelastic 
deformation (presented by inelastic strain tensor £'™). Therefore, strain energy 
f{E'^^) is obtained by decomposition of the function F{S,E*°*) into the "Taylor 
series" by the parameter E'^K This leads to the expression of elastic moduli of a 
media through the invariants of material metric g and the lapse function S. 
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6. Action, boundary term, Hooke's law. 

The Action functional is the integral of the Lagrangian density £(G, (f)) over a 
4D domain U = [0, To] x V. Here (V, dV) is an arbitrary subdomain of B with the 
boundary dV, combined with the 3D boundary integral that accounts for the work 
W of surface traction (0)- 

Au{,G,cl))^ f {L^{G)+Le{E-\G))dV+ f Vr{<jy,G)S^. (6.1) 

JU "'[0,To]xaVi 

Here d^S is the area element on the 3D boundary dU of the cylinder t/ (|S]). 

The second term on the right represents a boundary conditions put on the de- 
formation history 0. Typically the boundary dV of the domain V^ C -B is divided 
into two parts V = dVi U 91^. The deformation is prescribed on the part 9V2: 
4>dV2 — 4'itj^)j while along the part dUi of the boundary the traction r is pre- 
scribed. Function Vr depends (conventionally) on the velocity V of deformation 
and on the traction 1-form t, which is chosen in such a way as to have — V^Kr = ''' 

- traction (jEl)- In Euclidean space with the dead load one takes Vr = —t ■ (j). 

Deformation (/)(0, X) and the velocity V — 4>*i-§f of material points is assumed to 
be given at the moment t = 0. This determines initial conditions for the deformation 
history. 

The boundary conditions for the metric G (including initial conditions for 3D 
material metric g) require some special attention. Initial values oi S, g, N are known 

- prescribed by the material manufacturing process and by the previous history of 
the material deformation. On the part dVi of lateral surface dV where deformation 
4>\dVi is prescribed (for instance when this part of surface is not moving at all, 
see |19j for examples) we can find 4'*{g\dVi) by measuring distances between the 
material points on the boundary of the body in the physical space at moment t and 
recalculating them back to B by the tangent to the (prescribed) mapping 0( . If a 
part of surface is free from load, one can use the natural (Neumann type) condition 
that the mean curvature (with respect to the metric induced by gt) of this part of 
surface is zero. Along the part V2 of the surface where the load r is applied one may 
use for g analog of Laplace- Young condition for liquid surfaces relating difference 
of pressure with the surface tension and the mean curvature. The formulation of 
corresponding boundary conditions are the subject of another work. 

From the requirement that the variation of the action near the lateral sides of 
cylinder U are zero we get to the natural boundary condition 

P-N = T, or Y^aijPjN'' =T, 
1.1 

in terms of the first Piola-Kirchoff stress tensor P^ defined by the equation (material 
form of the Hooke's law, see [0]): 



dL, ^ df 






Notice that if the kinetic energy is included into L^ , Piola-Kirchoff Tensor has the 
density of linear momentum vector as its P° components (' PUll^ 'l. 

We will be using the second (material) Piola-Kirchoff tensor Sj = P^cp^j- 
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It is useful to recall that the (laboratory) Cauchy stress tensor aij is related to 
the first Piola-Kirchoff tensor by the following formula a-ij — J~^{(t>)hik(t>\Pj , J{4>) 
being the Jacobian of the deformation (f>. 

The zero condition for the variation at the top (T = Tg) a-^d the bottom (T = 0) 
of the cylinder lead to the relation between the linear momentum and the kinetic 
energy in the classical case (Legendre transformation). In the scheme presented 
here these variations also includes terms related to the aging processes. 



7. Euler-Lagrange Equations. 

The variation principle of the extreme action 5 A = taken with respect to the 
dynamic variables (j) and G results in a system of Euler-Lagrange equations that 
represent the coupled Elasticity and "Aging" equations 

^(Po«) + ^-Ea|7(^)-^"/N(VB)™ = 0, m = l,2,3. (7.1) 

'^^" ^^' /\G\T'\ I, J = 0,1,2,3. (7.2) 



SGij 6Gij 

The Elasticity Equations (7.1) are obtained by taking the variation SA with 
respect to the components (p^ within the domain U. In the case of a BD metric 
G and the synchronized deformation </>, these equations coincide with the conven- 
tional dynamical equations of Elasticity Theory. However their special features are 
associated with the different form of the elastic strain tensor E'^'' and with the de- 
pendence of the elastic parameters on time through the invariants of the metric G. 
The evolution of these parameters is defined by the equations (7.2) (referred to as 
Aging equations). 

The Aging Equations (7.2) resulting from the variation of action with respect 
to the metric tensor G describe the evolution of the material metric G for a given 
initial and boundary conditions. 

The right side of the equations (7.2) represents the (symmetrical) "Canonical 
Energy-Momentum Tensor" y/\G\T^'-' = ^ sg" ( ^ )■ -'^-'^ '-'^'' situation this 
tensor is closely related to the Eshelby EM Tensor bjj. 

In his celebrated works J. Eshelby ( \2'6l I15|). introduced the 3D and then 4D 
dynamical energy-momentum tensor (Eshelby EM Tensor) b (denoted Pj* in 

m)- 

b'j = fS'j E If 0:, = fS'j - 55, (7.3) 

f being the elastic energy per unit volume. 

The tensor b includes the 3D-Eshelby stress tensor f |15[l21| ') , the 1-form of quasi- 
momentum (pseudomomentum) V = 6j, J = 1,2,3, (see |15l I20p . strain energy 
density 6q = —Le = f (plus kinetic energy, if the last one is present) and the 
energy flow vector s ^ b^ ^ — P/^'q, / == 1, 2, 3. In the quasi-static case 6j = for 

J = 1, 2, 3. In the case of a BD metric (N = 0) G we have Vb =boB = 0, B ^ 1, 2, 3. 
Tensor bij is, in general, not symmetric (although its 3x3 space part is symmetric 
with respect to the Cauchy metric C3 (</>), see [201 ■ 
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It was proved in {22, that if metric G is block diagonal (i.e. if A^ = 0) and the 
body forces are zero, then 

where h^^'-''^ is the symmetrical part of the 4D Eshelby tensor and the symbol exp 
refers to the derivative of ie by the explicit dependence of G (not through S^'). 

For the Lagrangian L = Lm + Le defined by (5.4-5) the Aging Equations (7.2) 
can be rewritten in the more convenient ADM notations. 

The aging equations (7.2) take the form of the system of PDE for the lapse 
function S, shift vector field N and the 3D material metric g. The explicit form of 
the above equations can be readily obtained for the Lagrangian in a form (5.4-5). 
In order to achieve this the variational derivatives of components of Lagrangian 
with respect to the variables S, N, gu need to be calculated. In Appendix B we 
calculate the variations of some of these terms and present them in tabular form. 

Variation by S (assuming that / does not depend on S) : 

^ ^O^iF + S^) + ixiK)-^: K) + adiv,{Nf+m{9) - f + S^. (7.5) 
Variation by N^ : 



^ ^ 'iPP)"' - ^^"' ^ '""'°'> ■ ""''<"» ^ 1-''" J^J""' 



1 ^ /„ /TTT ^X AS^ ^ , 1 a t. fTZ\dx 



«svw'^'''°^''"3J^' '"' " ^^^^'^''^"^''"s?;" ^ »■ <""' 



Variation by gu (for simplicity, we omit in this equation the terms coming from 
diVg{NY in Lagrangian (5.4-5), for the corresponding term in the equation see 
Appendix B): 

Pa J 



aA7*"''"i«7-st* ' ;«^-s?j-' '^■'-;wki''r^''''^'' 



1 n f ^ /T-l 9X fA.rA, 1. AB . 9F 



rdx'< (/^ov^^s'^^'')] + ^^™5 



PO^/M V ^^B / 2 '"" dgAB 



d\\N\\l 2" "3- > 2'' "^ 2 dgAB 2 OgAB 2 

(7.7) 

Here S"^^ = Ric{g)^^ 2 9^^ ^^ ^^^ Einstein tensor of metric g. In the 

first line of equations (7.6) (left side) Ag is the 3D Laplace operator is defined by 
the metric g, Hess{f) = f:M;N stands for the Hessian of the function / (double 
covariant derivative tensor of /) . 

On the right side of (7.7) remains the symmetrized Second Piola-Kirchoff Stress 
Tensor S (Here and thereof S(^ab) = ^{Sab + Sba)) or the Eshelby stress tensor 
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since ^^^^^ = -b'^^^^ + Ceg^^ ■ The Eshelby EM Tensor b is thus the driving force 
of the evolution of material metric g (comp. |24jV 

Equations (7.1-2) together with the equation (3.4) for the reference density form 
a closed system of equations for dynamic variables (G/j,0*,po)- Complemented 
with the initial and boundary conditions, these equations provide a closed non- 
linear boundary value problem for the deformation of solid and evolution of the 
material properties. 

In general, the system (7.1-2) seems rather complex, especially if Le depends 
on the metric G an its (differential) invariants explicitly. Nevertheless, leaving a 
detailed analysis of this system for future studies, we make some brief remarks 
about special cases where system (7.2) is effectively simplified. 

8. Special cases and examples. 

8.1. Block-diagonal metric G. In a case of a BD-metric, iV = (no shift). 
Therefore, the metric Lagrangian has the form L,n = F{S,E'^"') + x[K) + (iR[g) 
that includes time derivatives of 3D metric g (in x{K)) and the space derivatives 
of g in the curvature term R{gt)- 

No derivatives of the lapse function S appear anywhere in Lagrangian. In par- 
ticular, equation obtained by variation of S is not a dynamical equation but 
rather a constraint, similar to the "energy constraint" in the Einstein equations 

(El)- 

In the case, when the elastic coefficients do not depend on S, equation (7.5) takes 
the form 

{F + S^) + (xiK) -§^--K)+ mig) = f, (8.1) 

where pof\/\g\S is the density of strain energy (per unit of unperturbed volume). 

This relation represents an equilibrium between the strain energy in the material 
(residual stresses presented by the scalar curvature of g) and the internal material 
constituents (the "ground state term" and the terms defined by the kinetic of 
material processes). In the case of a homogeneous tensile rod (^Hl) this relation 
determines the domain of admissible evolution in the phase space and the " stopping 
surface" where evolution of the material under the fixed conditions stops (see Sec. 11 
below) . 

As / ^ and the kinetic processes are stopped, the system tends to the "nat- 
ural" limit state which determines the relation between the "ground state energy" 
F{S,E'^'^) and the residual stresses (see Sec. 8. 3 below). 

8.2. Spacial subsystem. The spatial part (7.7) of aging equations represents the 
system of PDE for the metric gu having the form 

-/3f ^^(g) - S-'S'^Qfi'.ie.gVJ + W^ 5(^^). (8.2) 

Here Qf^f = g^^i qj^m and ^p — dt — N is the principal part of the 1st order 
linear operator ^ = dt — Cj^. The term W in the left side depends on the metric 
coefficients, function S and their first derivatives. Einstein tensor £[g) is linear by 
the second-order space derivatives of g. Thus, this system is quasilinear evolutional 
second order system for metric g. It can be easily transformed to the normal form 
under simple conditions on the dissipative potential x- 
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8.3. Statical case. Consider the case where N = Q,U = 0, f docs not depend on 
G,S exphcitly, S,g are time-independent, (3 = const. Then the system of aging 
equations is reduced to the following form (here and below F = SF) 

V^y^ - P^^"" + {pnS)-^[^,{p,S)g^^ + HessfB(poS)] + \{F + /3i?(g)).9^^ = i(/.9^^ - S^^^)) 

(8.3) 
In the absence of the strain energy, i.e. when f{E'^'') = 0, fe^'^-^' = system (8.3) 
has the trivial solution S = const, g = go- 

Calculate S^^ through the Cauchy stress tensor using (6.3) as follows: S^^ = 

kAaJ „CB _ ( T{A\^-i Auis„,\j,] „CB _ T(A\r.BCA.O w^-1 Auis^,, _ ^/\K^BC„A 



Multiplying the first equation in (8.3) by ^g"^^ and subtracting from the second 
we get 

' ^^ ' ^^g^^]^\ + {p.S)-'^\[/^,{poS)g^^+Hessj^{poS)]-P£^^^\ = -i/Rg^^'^a^^ 



'dgAB 2 83" JVI^MVru^; V I^IL-avru^;^ .^^^^^^g yrv^ n .-^^ Vl^l ^ 

(8.4) 
This is the balance equation between the metric characteristics (Einstein tensor, 
" ground state energy" , lapse function S) and the stresses in the body. It is especially 
simple in the case where 5 = 1 is absent from F: 

^^ /TTT RpAB,/r-\, ^-1. /Or/A ^.\„AB . TT^..AB„, _ _l . /Tu JB\C A) 



l3£^^^\+p,'^\[{Agp,)g^^ + Hessfpo] = - ^\g(B\^a', 



^ V i.yi A^*-- V lyi I Po V i.yiiv^spivi' 1 -1-11.^.3^ puj — ^v I'^y ^c ■ 

(8.5) 
Here we can see how the curvature of material metric and the density of non- 
homogeneities may be a source of the stresses in the body in the absence of elastic 
deformation, i.e. when the conventional strain tensor E^^ '^°" = ^ln{gQ^C3{(j))) is 
zero. Namely, in such a case though the conventional strain tensor is zero, decline 
of the Cauchy metric C^{(t)) from the material metric g is not zero. Subsequently 
stress tensor S is not zero. Equation (8.5) thus describes the self-equilibrated stress 
resulting from the curvature of the metric g and is related to the incompatibility of 
embedding of the solid into the physical space. The first term on the left in (8.5) 
is related to the deviation of the total energy from its stationary value. 

One example of this situation a nonhomogeneous chemical transformation (ox- 
idation) of material, which results in the variation of material density and an in- 
compatibility with the reference configuration. A more specific example of stress 
induced chemical transformation is discussed below in section 8.6. 

8.4. Almost flat case. Here we use essentially that the dimension of B is 3. In the 
case, where Ric{gt) « 0, a good approximation of the general system (7.1-2) can be 
proposed. If the total deformation cj) is approximated by the " ground deformation^' 
<j) (i.e. deformation (j>{X,T) such that (j)*h = gr, recall that this is the synchronous 
case!) in the evaluation of the EMT T^'^ on the right side of aging equations (7.2), 
the latter becomes decoupled from the equilibrium equations (7.1). This allows 
us to study the aging equations separately from the elasticity equations and, after 
obtaining solution for G, substitute them into the elastic equilibrium equation (7.1) 
and solve it as the conventional elasticity equation with variable elastic moduli. 
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8.5. Homogeneous media. In the case of a homogeneous material (|20]) metric 
G depends on T only, and Einstein tensor £ij{g) is identically zero. As a result, 
(7.2) becomes a system of quasi-linear ordinary differential equations of the second 
order for the lapse function N and the material 3D metric gij. The Cauchy problem 
for this system is correct under some mild conditions to the dissipative potential x- 

The linearized version of aging equations of ID homogeneous rod was discussed 
in ([2S])- In Sec. 11 we shall briefly describe the study of some aging problems for a 
homogeneous rod (more detailed presentation will be published elsewhere, see (19p. 

We conclude this section with two model examples that show the type of material 
behavior that can be studied using presented approach. 

Example 1. Modeling of Necking Phenomena in Polymers. 

Delayed Necking, observed in various engineering thermoplastics, is a pictorial 
illustration of traveling wave solution. Necking in general is a localized large defor- 
mation (drawing) of a polymer with a distinct boundary between the drawn and 
undrawn material domains ( |26II27[I^ ). Delayed necking takes place in a rod in 
uniaxial tension, i.e., under constant applied load when the initial Piola-Kirchoff 
stress 511(0) is less then the yield stress. At first a uniform creep takes place, 
i.e., a uniform stretching with a draw ratio A — I/Iq, where I stands for an actual 
(current) length scale. After certain time interval when the increasing stress 6*11 
reaches the yield stress value, necking, also called "cold drawing" with a natural 
draw ratio A = h/lo starts, i.e., strain localization is formed and propagates along 
the rod with a constant speed N^. The observed elongation results exclusively from 
a transformation of the original material adjacent to the neck boundary into the 
drawn (oriented) state and propagation of the boundary along the rod, as depicted 
in Fig. 6. 

We consider here a ID model of a rod, with the lapse function S = I and the 
shift vector field N = N^ -^ being constant (see JZj for a 3D model of the necking 
process). Denote hy g = gii(t,X) — X^go the only component of material metric. 
Take the "ground state" energy as 

F(A) = (A-Ao)2(a + 6(A-Ai)2) 

where the elongation of the rod X{t, X) = i— is the drawing variable, Aq — 1, Ai 
are two states (to compare with example of the creep in Sec. 1 1 put A = e** w 1 + r]) . 
This is the simplest function that admits two different stable states (metrics) with 
equal chances when true stress reaches a critical value. 

The metric Lagrangian is reduced to Lm = F{\)+x{K), where K — {g~^uc-g) — 
29t-77, 77 = ln{\) where dr — uq — dt — N^dx- 

Experimental data suggest that in the necking the material density variation is 
negligible, thus we take pQ{t,X)\{t,X) = pq{0,X). As a result, the action takes 
the form 

A{\{t,X))^ f [F{X)+xi'2drln{\))]dtAdX = 

J[o,ti]x[a,L] 



[(A - Ao)'(a + h(\ - \xY) - —(A - Ao)" + x{'^drln{\))\dt A dX, 

[0,ti]x[0,L] ^0 

(8.6) 

where ^0 is the initial cross section of the rod and J- is the force acting on the right 
end pulling in X-direction. The second term here represents the work of the load 
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Figure 6. Necking of a ID Rod 

on non-elastic deformation. Consider the case where J^ is large enough to change 
sign of the quadratic part of the" ground energy" F. For simplicity we take -j- = a. 
Variation by A leads to the aging equation in the form 

xKr - yk' + (2A-iA.)x"] - ^^'(A) = 0, 

where F(A) = h{\ - \of{\ ~ \if . 

The 2D dynamical system corresponding to this equation has equilibria points 
(Ai,0), « = 0,1,2 at the roots Ao,A2 = ^(Ao + Ai),Ai) of the polynomial F'{\). 
If the dissipative potential x(m) satisfies to the conditions x'(0) = 0,x"(0) > 0, 
root A2 is the center while other two are saddles whose separatrix loop enclose the 
elliptic region. 

For a given F{g), when the stress reaches the initiation level and from the trivial 
solution for g there bifurcates the separatrix solution, then, we get the "traveling 
wave solutions" in the form of "kink" (^7]), propagating with the the speed N^ 
along the rod, for the metric g{X,t) (and, by the mass conservation law, for the 
density po(-'^,0) ■ 

8.6. Example: Variation of Material Metric g due to the Chemical Degra- 
dation. Here the evolution of the uniform material metric to the piecewise constant 
metric with the jump along the interface between a layer of chemically degraded 
material and the original material follows the kinetics of chemical degradation (see 

my 

Consider a thin- walled thermoplastic tubing employed for transport of chemically 
aggressive fluid. In time, the inner surface layer of material undergoes chemical 
degradation due to interaction with aggressive fluid flow. 

Chemical degradation is manifested in an increase of the material density po, 
significant reduction in toughness (resistance to cracking) and a subtle change in 
yield strength. Young's modulus and other thermo-mechanical properties. 

Assuming the homogeneity of degraded layer we see that the original euclidian 
material reference metric in degraded ring evolves (see the mass conservation law 
(3.4)) which generates a jump on the interface with the outer layer of unchanged 
material. Continuity of normal stresses on the interface allows us to describe the 
final state of the system by elementary methods presented below. 

Consider a thin ring (see Figure 7) which represent the 2D cross-section of the 
tubing. The wall thickness t — Ro ~ Ri is small in comparison to the outer radius 



22 ALEXANDER, CHUDNOVSKY AND SERGE PRESTON 

Ro'- t/Ro <C 1. Rd in Fig ** stands for the radius of interface between the layer of 
degraded material and unchanged layer. The depth of degradation td = Rd — Ri is 
relatively small: td/t ^ 1. 

Select the polar coordinate system (r, 6). 2D material metrics of the initial {g^) 
and degraded {g') states are 

9' = {lr(.f)AA-r\Q);g'^{'l';^)A9'\ = {l + e)r\ 

where r' = (1 + e)r and e is a small variation of scale in the radial direction. 

Mass conservation law poylff^ — p'ayW] relates density variation Pg = po + 
ApOj —^ ^ 10^'^ with the change in material metric 



Po 



- = J^^^il + e)y^=^l-^^l + le + Oie^). (8.7) 
'0 y r(Oy po 2 



Po + Api 

Therefore e w -|^^ + Ofe^). 

Thus, the densification (i.e. Apo > 0) leads to the shrinkage of the thin ring of 
degraded material. If we remove the constrains on shrinkage applied by the outer 
ring of original material the gap 

^-R"d-{i-l^)R'd = l^R'd (8.8) 

3 Po J Po 

appears. 

As a result of such constrains, the degraded material should be elastically stretched 
to close the gap w. This elastic deformation has the form 

ir"(r) = (1 + ^^i^)r' 



Under such a deformation elastic strain tensor E'^^ = ^g ^(C3(0) — g),Cz{<t>) 
hij4)^ jcfp J has the form 

e/., 2Apo A 



The tensile strain (8.9) is directly translated into the tensile radial stress via 
Hooke's law 

Y 2 Apo 

1-1/3 Po 
Although hoop stresses are may be discontinuous, the equilibrium conditions re- 
quires continuity of radial stress across the interface, , 

Crrnrlr^Rd-O — "'rr^r |r=i?d+0- 

This implies that the outer ring of original material experiences compressive stresses 
while the inner degraded layer is under tension. The elastic strains i?*^' jointly close 
the gap w and restore the compatibility of the Cauchy metric g final in the whole 
domain: 

g/znai = g' + E''K 

Therefore while the material metric g' has the jump leading to the nonzero singular 
curvature along the interface, the final metric is continuous and flat. 
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Figure 7. A sketch of polymer tubing cross section with inner 
degraded layer 

9. Physical and Material Balance Laws 

As it is typical for a Lagrangian Field Theory, action of any one-parameter 
group of transformations of the space P x M, commuting with the projection to 
P, leads to the corresponding balance law (See j^])- In particular, translations in 
the "physical" space-time M lead to the dynamical equations (7.1-2), rotations in 
M lead to the angular momentum balance law (conservation law in the absence of 
applied torque). Respectively, translations in the "material space-time" P lead to 
the energy balance law (translations along the time T axis) and to the material 
momentum balance law (" pseudomomentum" balance, ([211, POji IH^); rota- 
tions in the material space B lead to the "material angular momentum" balance 
law ([2n|)- 

In the table below we present basic balance laws together with the transforma- 
tions generating them. It is instructive to compare the space and material balance 
laws as it has been considered previously by several authors ('|3(Jj. |31j). 

Table 1. Space and Material Balance Laws 



Physical space-time 
(Material independent) 



Material space-time 
(Space independent) 



Symmetry 



Linear momentum balance law 
(Equilibrium equations) 

divjcr) = f 



Material momentum (pseudo- 
(-momentum) balance law 
div(b)=fmat 



Homogeneity of 3D-space 



Energy balance law: 



Time homogeneity 



Energy balance law 
iVlaterial angular momentum 
balance law = C-symmetry of 
Eshelby stress tensor b 
h:C = C:h 



Isotropy of 3D-space 



Angular momentum balance law 
= h-symmetry of Cauchy stress 
tensor a: I : cr = cr : I 



Space and Material balance (conservation) laws are related via the deformation 
gradient d(j). Restricting ourselves to the synchronized case and writing the material 
balance laws in the form rji — and their " physical" counterparts in the form Vi = 
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we get the relationsship between these famihes of balance laws 



Vo 



\ ■■■ <^\ 



1 ' ,. 

A3 



\ (9.1) 



V ct^% J 

Similar to the case of relativistic elasticity (J3), the system of material balance laws 
rji = 0, 1—1,2,3 is equivalent to the elasticity equations i^i — 0, i — I, 2, 3, while 
the energy balance law rjo = (which here is the material conservation law as well 
as the physical one: in the case of synchronized history of deformation (f> material 
time T and physical time t coincide). As a result, the energy conservation law is the 
consequence of the time translation invariance in both senses and follows from any 
of these two systems: 770 — J^lZi 4'^o^i- This reflects the fact that the deformation 
we consider here are not truly 4-dimensional. 

Balance laws (with the source terms) can be transformed into conservation laws 
by adding new dynamical variables. In the theory of uniform materials ( 20, 24, 32j) 
it is zero curvature connection in the frame bundle over M that is added to the list 
of conventional dynamical variables, in our scheme - it is the 4D material metric G. 

10. Energy-Momentum Balance Law and the Eshelby Tensor. 

In this section we consider the Energy-Momentum balance law resulting from 
the Least Action Principle and the space-time symmetries. 

Consider local rigid translations in the material space-time X'^ \ — > X'^ + SX"^ . 
They generate a variation of components 0' of the deformations, components Gjj of 
material metric and their derivatives (we follow the arguments of J. Eshelby (',21,). 

Taking the variation of the Lagrangian density C — ^/\G\L — ^/\G\{Lg{G) + 
Le{G,E'^'-)) with respect to the material coordinates X'^ , one obtains 



i—0 i—1 



dx7 r ^^^'^ - dG^^^-''' + dGJ^^^^^' " dX^ 1 9G4f ) ^^'j ^^°-^^ 



The last term in the right side of (10.1) includes a definition of the (l,l)-tensor 
density £{G). Employing the Euler-Lagrange equations (7.1-2), we obtain for the 
Total Energy-Momentum Tensor (density) 

£*°* = -£^5_S'5+5(G)S, (10.2) 

the conservation law 

dlVGo (f *"*) = gJ^^T ^ - 0, J - 0, 1, 2, 3. (10.3) 

Divergence here is taken with respect to the 4D "reference" metric Gq. Since L^ 
does not depend on deformation (j) and the body forces potential U does not depend 
on its derivatives while Lg = —pof — PoU, 

'55 = - E S&€j = E P'^JS^/Wl (10.4) 

which is the 4D-version of the (density of) Second Piola-Kirchoff Stress Tensor. 
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Rewrite f *°* in the form: 5*°* = B + U5^jy/\^\+Cm6^j+E{G) where we denoted 
by B the tensor density B = b^J\G\ of the Eshelby EM Tensor. Then the equation 
(10.3) takes the form 

dlVG.B = B'jj = -diVGo{US'j^/M\+^rnS'j+£{G)), (10.5) 

where in the right side only metrical quantities and the potential U of the body 
forces are left. 

The second term and the metrical part of the third term in the right side of (10.5) 
are related to the ground state of the Lagrangian density i.e. to the inhomogeneity 
of "cohesive energy" and the "material flows". The elastic part of the third term 
on the right is related to a variation of elastic moduli if these moduli depend on 
the derivatives of the metric G. Equality (10.5) can be easily rewritten in terms of 
covariant derivatives with respect to the metric G. 

Taking J = in (10.6) we arrive at the energy conservation law in ADM notations 
(using t instead of X^) 

d ,,, ,. . _ ^ „n. 4^ a ^'=^ ^ 



((/ + C/ + L„)5/M+0 = E77w E^''^:o^^-^(G)o- (10.6) 



Q^KKJ ' - ' -™/-V,.., ' -0/ A^QXI 



= 1 



Equation (10.6) has the form — — 2_2 — 2^1 — "^"^^ y> — TotalFlowDensity, with 
the total (inner) energy density given by 



£'o'"° = if + U + L^.)SV\9\+£iGfo- (10.7) 

The total energy is the sum of the following parts: elastic energy /, potential energy 
of the volume forces U , cohesive "ground state" energy - the term F{E'^'^, S) in L,„, 
inhomogeneities energy from the curvature density f}R{g)S^/\g\ and corresponding 
terms of £q , "kinetic metric energy" that is defined by the term s produced by x{^) 
in L,n and £q and reflects the intensity of irreversible deformation and "metrical 
volume change energy" coming from the diVg(iV)-terms. 

The sum on the right side of (10.5) consists of the flow of the Piola-Kirchoff 
stress tensor density Y^jZi[Pi4'''o'^\A9\]j ^^"^ H^^ flows related to the change of 
the material metric - internal material flows, flows of inhomogeneities (coming from 
the curvature R{g) etc. 

If the metric G does not depend on time (i.e. N — 0,K — 0, G = Go) and 
if Ric{gt) = 0, one obtains the conventional energy conservation law of Elasticity 
Theory (Ej, Chapter 5, Sec.5): ^^ = -E'iZI o^iP^lo),- 

Example 2 (Block diagonal metric G, synchronous deformation and homogeneous 
media). In this case we have Ric{gt) =0, iV = 0,.g = g{t), S = S{t), the extrinsic 

curvature has the form K{ = \ ^ n ^ tk I and \(K) is the only term in 

the Lagrangian containing time derivatives. 

In addition to this, no flow terms except the usual Piola-Kirchoff flow appear on 
the right side in the energy balance law which takes the form 

^ ((/ + U + L^)sVW\ + S-'gABM [1^/'^ + ^a"""] SVW\^ =+J2(P'^:oSV\9\) 

'"(10.8) 



dt 
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This equation describes how the energy supphed by the boundary load spreads not 
just to the increase of the strain energy, but also to the change of its "cohesive 
energy" of the material (^^(5*, \g\)) and to the acceleration of the aging processes. 

11. Aging of a homogeneous rod 

In considering three types of inelastic processes in a tensile homogeneous rod: 
unconstrained aging, stress relaxation and creep (see |19| for more detailed exposi- 
tion) we assume that N = 0, R{g) = 0, S{t = 0) = 1 and that S{t) is increasing to 
a certain level depending on the initial state of the body and the process. 

11.1. Deformation, strain tensors and tensor K. Introduce material cylin- 
drical coordinates {R, Q, Z) in the reference state of a rod B. Spacial cylindrical 
coordinates {r,d,z) are introduced in the physical space R'^. In addition we nor- 
malize the initial state g(0) of material metric taking g{0) — go- 

We consider the class of time dependent (total) deformations (j) of the from 

0t : {R, e, Z)^{r = ^j.{t, Z)R, 9 = e,Z = k{Z, t)), (11.1) 

with fi, k representing amount of "stretch" in radial and axial directions respectively. 
Material metric gt is flat (homogeneous case!) and is generated by a global 
deformation 0™ of the same type as (11.1), with firm ^m the same as above: gt = 
(f>jn{t, -yh. As a result 

g = €ah = R'A. , C(0) = 

■(11.2) 
where A — kz^^m = km,z- For a homogeneous rod k{Z,t) = X{t)Z, km{Z,t) — 

Xmit)Z. 

In the elasticity theory (see |22|) it is customary to present deformation (total 
and inelastic as well) as the composition of a uniform dilatation with the axial 
expansion factor A„(i) and of the volume preserving normal expansion with the 
factor A^: \ — XyXd- We will obtain fi ~ XyX^ and \/\g\ = X^R (and the same 

for Xmv, XmdJ- 

As a result, the inelastic strain tensor can be written in the following form 

1 , /in(p) \ fi-hn \ 

2 V hi{\}J V e+r,/ 

Here we introduced the variables ^ = ln(Xm v),'n — ln{Xm d)- As a result, 
calculating basic invariants of these tensors we see that the " ground state" energy 
F is the function of 5, ^, 77^. 

Decomposing the total deformation as the composition of inelastic and elastic 
one and assuming that elastic deformation is small compared to 1 we write: 

At, = At, m(l -f tv), Xd — Xd m(l + Cd)- 

In these notations elastic strain tensor takes the conventional diagonal form 

E"^ = -jj-nig-^Cicj)^"^)) « diag{ e„ - -e^, e^, - -erf, e^ + erf). 
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Strain energy for our (homogeneous) rod will now take the form 

f = §4 + ^el (11.4) 

with the bulk coefficient K and the Lame coefficient /i. 

Mass conservation law (3.4) takes the form po{t) = A^po(O). 

The spacial part of the tensor K for a homogeneous rod takes the diagonal form 

K = S-^diag{2^t - Vu 2^ + 277*, 2^* " Vt)- (11-5) 

Thus, Tr{K) = 65-1^*, Tr{K - \Tr{K)lf = eS'^r]^ and dissipative potential 
xiK) is the function of arguments S^^^t^ S^^rjt. 

Remark 4. In general, aging equation for the described situation have the form of 
a 3D degenerate Lagrangian system (we refer to ^H), or Pl] for more details). In 
the cases of the processes studied below this system reduces to the 2D degenerate 
dynamical system. In all three cases one can trivially solve elasticity equations, 
exclude elastic variables ev,^d from aging equations and, therefore, to close the 
system of aging equations. 

11.2. Unconstrained aging. Unconstrained aging (shortly UA) is the simplest 
example of a material evolution. A sample of material (rod) is prepared and then 
is left without any constraints or load applied to it. Usually the process of aging is 
manifested in a variation of material density, or a specific volume change up to a 
saturation point, when the observable evolution stops. In many polymers the aging 
is accompanied by shrinkage up to a few percent of initial volume. This diminishing 
in volume (2-5%) is called unconstrained aging. We discuss here a model for UA 
in terms of variables (5,^) (dilatational deformation plays negligible role in UA). 
No strain energy is present, stress is zero. 
We take the "ground state energy" to be 

Fua{S,0 = (ci + C2S + {p^S-^ + ke) 

with fc > 0,p < 0, ci < 0, C2 > and the dissipative potential 

Integrating over the volume of the rod we get the action in the form 

A(e, S) = V j {{ciS + C2S^ + (K + k^^S) + aSiS-^^tf)dt. (11.6) 

Jo 
Euler-Lagrange Equations of UA can be reduced to the following dynamical system 

~ " ^ ^ : (11.7) 

* 2C2 



We have here $* ^ 0, S"* ^ 0. 

Equation (7.5) takes here the form a^^ = S^-ka^, where F = SF. 
Take a = — 1. Then, the domain of admissible dynamics defined by the positivity 
of expression under the square root is 

ss-^{ci+ke), (11.8) 

2C2 
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Figure 8. Aging curves 

and the curve where evolution stops when the phase trajectory reaches the final 
state is S" = ~^{ci + fc^^). 

The " ground state energy" F is negative at initial moment and that it increases 
during the evolution. 

System (11.7) has the first integral J = 025*^ — p^. Choosing an initial point 
(S'(0),^(0) = 0) of a trajectory F in the domain of admissible dynamics. Along F 
we have J{£,,S) = J{S{0),0) If we calculate S{t) as the function of ^(i) along F, 
substitute into the first equation (11.7) and separate variables in this equation we 
get the ^(i) as the explicit function of parameters of the problem and initial value 
S{0) in terms of elliptic functions (see jJH). 

Figure 8 shows a family of shrinkage curves corresponding to the various values of 
S'(O) = 1, 1.3, 1.6, 1.9 (which represent the initial aging of the material). Apparently, 
the higher is the initial age, the less shrinkage is observed. 



When a load applied to the rod reaches certain level, new processes may start. 
These new processes (going on the background of the UA) initiate action of a new 
part of the "ground state" F{S, £,, rj) and activates the new kinetic potential X2{K). 
For the description of stress relaxation and creep we choose the dissipative potential 
X2{K) = -Mjin{%) — ^(1 + j))^n{l + ^) corresponding to the phenomenological 
Dorn relation between the stress and the strain rate rj (see [35], Sec. 2. 3 and the 
footnote in Sec. 5 above). Unconstrained aging is much slower and leads to smaller 
changes then both stress relaxation and the creep. That is why we may with good 
accuracy disregard the UA while describing two other processes. 

11.3. Stress relzixation. In the case of a stress relsLxation (SR) we fix the rod 
of initial length L at the left end and then quickly pull (or compress) it uniaxially 
and fast (elastically) until it reach certain length L*. Then we fix right end as well, 
leaving side surface of the rod free. In this configuration the only component of 
Cauchy stress that is nonzero is cTzz- For the SR the volume change is negligible 
and we have A,„ „ = 1 , A,„ = \„i d 

Initially all the stretching is due to elastic process and AJ = L* / L = (1 + 
ez(0)). Then the inelastic deformation starts to increase in expense of elastic one 
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maintaining the total strain constant. The reduction of elastic strain is directly 
translated into the reduction of stresses via Hooka's law. The total elongation at 
moment t can be decomposed as follows 

AS - (1 + e,{t))XrnAt) - (1 + ez{t))e^^''> (11.9) 

and therefore ez(i) = ^7* — v{'t)> where ij* — ln{X*). 

From Hooke's law azz ~ Ye^, where Y is the Young module. Thus for the strain 
energy expression we obtain 

fiv) = l<^zz^z = ^Yel = |(ry* - r^t))' ■ (H-IO) 

For pure stress relaxation (without background UA) 
F = FsrIS, r]) = (qi + q2S) + ry(6oS'"^ + &i + an]), qi < 0, 92 < 0, 60 < 0, ai < 0, 

with the coefficients different from those of the slow UA. 
Action A{S{t),ri{t)) now takes the form 



Jo 



\qi + q2S) + viboS-' +b,+ air?) + j{v*- vitjf + Sx{S-\) 



Sdt^ 



f \FsR{S,7j) + Sx{S-'vt)]dt, (ll.fl) 

Jo L J 



where 



Fsr{S,ii) = SFsR = {qiS + q2S^) + v{bo + biS + ai^S) + S^{,f - r]{t)f = 



q2S'+P2{ri)S, (11.12) 



where ^2(^7) = [qi + \v* ^) + (&i - Yjf)i] + (ai + ^)7f ■ 

We have Fs — 2q2S + P2{il) ^^^d the domain of admissible motion is defined by 

Da, = {{v,S)\v>0,S^l, S<-P^}, 

^92 

while the stoping curve has the form S ~ — ^^j^. 
Aging equations (7.5-7.7) will now take the form 

Ut = S^-\2q2S + P2{v)) = Si,-H2q2S + [(91 + ^77* 2) + (b^ _ y^*),, + (^^ + Y)^2^)^ 
\St = ^^-H2q2S + P2iv)) = S^-\2q2S + [{qi + \r,* ^) + {h - Yrj*)^ + (a, + ^)rj^]). 

(11.13) 
This system has the first integral J = boT] — q2S^, and the phase trajectory cor- 
responding to the initial value S = 5'(0), 77(0) = has the form S = Js{0)^ + ^ry. 

Using this one can find analytic solutions ri{t) in terms of elliptic functions ( fl9j ) . 
On the Figure 9 we present results of calculations of the stress relaxation azz (t) 
for several values of initial stretching L = 0.1, 0.15, 0.2, 0.5 and for realistic values 
of parameters of the problem. Values of azz{t) are found by solving numerically 
system (11.13) for r]{t), calculating elastic strain €z{t) and using the Hooke's law. 
Apparently, the higher is the value of initial stretching, the sharper is the stress 
relaxation curve and the higher is the asymptotic value of stress. 
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SR Figure 2 stress relxation 




Figure 9. Stress Relaxation: Stress azz{t). for t]* ^ 0.1,0.15,0.2,0.25. 



11.4. Creep. In the case of the creep we fix the left end of the rod and apply 
force T in Z-direction to its right end. If this force is large enough (i.e. if the 
concentration of elastic energy / in the rod is larger then an activation threshold), 
the creep starts - inelastic deformation that goes on for some time until the rod 
brakes. Thus, at the moment when the inelastic strain r](t) starts growing from 
zero, there should be a supply of strain energy obtained from the work of the 
stress cTzz 
initiation) 

During the creep the homogeneous component (Tzz of stress is equal 

T _ T _ \d{t)T _ e^T 



-^ on elastic deformation. Denote by /i„ this strain energy (of 



.w 



A{t) 



Ao 



Ao 



(11.14) 



The last equality is true given the assumption (natural for a conventional creep) 
that inelastic volume change is negligible, i.e. At, = 1, ^ = 0., for a constant 
force T and variable cross-section area A{t). 

Using the Hooke's law one can show that e^ = ^ and that the strain energy is 

equal to / = 2YA(oyi ■ 

Calculating the work of the load JF on the total way L{e^^^^ — 1) of the right 
end of the rod we get the additional term in the Lagrangian (work of load on the 
inelastic deformation) equal -j-ie-^ — 1). 

Overall action takes the form (^(i) — 0) 

C{S{t),T^{t)) = [Fcr{S,ii) + ^(e" - 1) + 4>{S-^rit) + ^T^e''^]Sdt, (11.15) 

where F{S, rf) is the same as for stress relaxation. 

Acting as before we get the following dynamical system for parameters 77, S: 



r,t = S^J-\^{e^ - 1) + AT^e^^ + (gi + 292^ + &i7? + a^rj^)), 



S, 



t = |7^"'(£(e'' - 1) + A^'e2'7 + (ft + 2^2^ + b,rj + a,7^^)). 



(11.16) 



where ip^x) 



^Dx 



1)+ function as in the case of stress relaxation. We have in 



this system r/t ^ 0, S** ^ for the admissible initial values 5(0), ?7(0). 
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For the creep to start, the argument of the function iIj{x) m the system should 



be positive at initial moment. Since (Tzz — ^—, this condition takes the form 



2Y 



qi+2q2S{0) + air]{0f>0. 



(11.17) 



Since qi, 52, oi are negative parameters, the inequality (11.17) defines the stress (or 
strain energy - /.„) threshold for the initiation of the creep processes (see |35, , p. 7). 
On the Figure 11 the graphs of ?7(t) for the creep are presented for different values of 
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Figure 10. Creep: inelastic strain r]{t) for F = 1, 1.5,2 and 3. 

the force -^ — 1, 1.5, 2 {Aq = 1). There is a point on each trajectory, corresponding 
to the instability of creep deformation where the cross-section of the rod diminishes 
to zero and the rod fails in so called ductile manner. As it can be seen from the 
results of calculations, the higher is the force, the faster creep deformation develops 
and the time to the ductile failure becomes significantly shorter, for example, three 
times increase in force results in more then 3 order of magnitude in time to failure. 
Comparing these graphs with the experimental data ( J35ll3f)j ') we see the good 
qualitative (and, for some materials, quantitative) agreement. 

12. Appendix A. Strain energy as a perturbation of the "ground 

STATE energy" . 

In this section we discuss perturbation scheme of a ground state Lagrangian 
F(£'™, iS*, TV), E™ — i \n{gQ^g) by elastic deformation, assuming that elastic strain 
tensor is small in compare to the inelastic one. 

In the pure inelastic mode of behavior (free aging, special load that produces 
C{(f>) = g) only metric quantities (7, S", N enter the total Lagrangian L. Under the 
general load the material metric g is deformed into C(0) and total deformation 
£:*°* = i \n{gQ^C{(j))) takes the place of E^™. 

Assuming that g~^C{4>) ~ I+small, so that E'^^ ^ £'™ we decompose gQ^C{(t>) = 
(go^g) ■ i9~^C{(j))). Taking logarithm, we obtain 



E*"* = - \n{exp{2E''' 



exp{2E''^)) ~ - ln{exp{2E'") ■ {I + 2E^^)) 
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Here we've used the linear approximation exp{2E^'') ~ I+2E^''. Using the Campbell- 
HausdorfF-Dynkin formula we get 

£;*"* ~ -[\n{exp{2E'"j) + Ad{exp{2E'"j)2E''^] = E'" + Ad{exp{2E'"))E''K (12.1) 

We consider the "total ground state" Lagrangian F depending on E*°* through 
its invariants /fc(i^*°*), k — 1,2,3. Then we decompose it into Taylor series by 
considering E'^'' small in compare to E™. The approximate expression for i5*°* 
above gives us up to the second order terms 

IkiE'"*) ~ Ik{E"') + dIk{E'''){E''^ 3^''i)+d^Ik{E"'){E''^ 3o'9^E^i ao'a^ + h.o.t.. 

(12.2) 
Substituting this into the metric term F{Ik[E*-°*-) (for this discussion we suppress 
in F all other arguments it depends on) we get, after recombining its terms, its 
decomposition up to the second order 

F(£:'°*) ~ i^(£;") + ^F/, (£:'") < dIk{E'"),E''^ 3o'9 > + 
fe 

+ {^FjkiE"') < d^ Ik{E"')[E''^ So'a^E"^ 9<7'9)+ 
fe 
+ ^F/,/^.(£:") < dIi{E"'),E''^ so'a >< dIj{E''^),E''^ So'9 >} (12.3) 

ij 

In this formula first term on the right is the basic metric (" ground state" ) energy 
describing, in particular, equilibrium values for the metric g (see example below). 

The second term, linear by E^^ describes the interaction of elastic and inelastic 
processes. In the absence of such interactions or other material processes g takes the 
value delivering minimum to the basic energy F{E''^). Therefore, its differential 
takes a zero value for corresponding value of the argument -B*" and linear term 
vanishes. 

Finally, the quadratic form in formula (12.3) is the conventional elastic (strain) 
energy with variable and possible inhomogeneous elasticity tensor. 

During the active processes of configurational changes in the material (aging, in 
the zone of phase transition) this elasticity tensor as well as the basic energy plays 
an active role in the evolution. But when such processes stops (no aging happens 
or wave of phase transition passed) and the material metric g is locked in some 
stable state (local minimum of F7), the value of this tensor is also locked at the 
corresponding value (see below) . 

In order to calculate the elasticity tensor e we have to calculate differentials of 
invariants Ik of the inelastic strain tensor i?*". We choose momenta Tr{E^) as 
the basic invariants of a (l,l)-tensors (|321)- In our, 3D case we have Ii{E) = 
Tr{E),l2{E) = Tr{E'^).h{E) = det{E). Thus, we have for its first diff'erentials 

said <dIi{E),P>=Tr{P); < d^iE), P >= 2Tr{EP); < dh{E),P >^ iTr{E^ P). 
Here we are using multiplication of (l,l)-tensors. The second differentials of mo- 
menta have the form d'^h{E) = 0; d'^l2{E){B,B) = 2Tr{B^); d'^h{E){B,B) = 
QTr{EB'^). 
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Before using these expressions for differentials we notice that since g^ g — 
exp{2E™) , (_E*")9o 9 — £;«" ^^lA due to the properties of Tr for all natural powers 
a,b one has rr(^*" °'[[E^i')9o'a f) ^ rr(^" "E''^ ^). Thus, conjugation by g^^g 
disappear from the formulas for Elastic energy. 

Substituting expression for differentials in (12.3) we get expression for F{E^°^) 
as the sum of " constant" , linear by E'^^ and quadratic by £''^' terms 

F(£:*°*) ~ F(£;™) + Tr{CE^^) + Tr{e : E"^ : E"^) (12.4) 

Here C is the (l,l)-tensor 

C = Fj,{E'")I + 2Fj,{E'")E"' + 3Fj,{E'"){E"')^, (12.5) 

and e the elasticity tensor 

efg = 2^,,,(£;")5g<5^+6F,3(£;™)g5^+2F,,,,(£;™)<5f5g+4F,,,,(i?")f(£;™)g+ 

9i^,/3/3(^")i''(^^")c''+4J^././.'5f(i?^")g+6F,,,3<5|(i?-)2,^+12F,,,3(S™)f(i?™)2,^. 

(12.6) 

If the ground state function F is given as the function of variables g, S, N, these 
formulas determine values of elastic moduli in a material which depend on the 
point X and on time t through the metric variables g, S, N. If these variables 
take stationary values, we get isotropic but nonhomogeneous material. If they are 
constant - we return to the conventional linear elasticity. 

Example 3. Isotropic material For isotropic material of linear elasticity, elastic 
tensor (in its (l,l)-version) has the form ( 6 ) 

eii = 2^^5l6i + X6i6l (12.7) 

Comparing with (15.7) we see immediately that there are two simple cases when 
(15.3) determines an isotropic material. 

Case 1 - generic. Take Eg^ '^ ~ hS^ where h is a. scalar function of g, S, N. 
In this case 

+ [2Fj,i, + AFj,i,h^ + QFj.iX + 4=Fj,i,h + 6Fj,i,h^ + 12Fi,i,h'']S^6S ■ 

(12.8) 

The first bracket gives expression for 2fi while the second one - for A. 

Case 2 - simple elasticity. In this case we have no aging, i?"' = 0. Then we 
get material with 

eiS = 2FjMSES^ + 2Fj,jMSiSS. (12.9) 

Thus, in this, restricted case 2/i = 2_F_/2(0),A = 2^/^/3(0). 

Example 4. Consider the model ID case with one component of strain tensors 
j^ei ^ j^in ^ ^j-jvial decomposition i5*°* = E'^'' + E™ and simple Taylor decomposition 
of the basic energy function F{E*°*): 

p^^tot ^ pi^Ein^ j^ F'{E'")E''' + -F"{E"'){E^^f + h.o.t. (12.10) 

As a result, strain energy here has the form 

[/(£;"') = -F"{E'''){E''^f, (12.11) 
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and the Young's module (or compressional stiffness, in a case of an elastic bar) is 

y = -F"{E'"). 

Consider two special cases. 

1. Classical elasticity. In this case we take F{E) = Fq + cE^. In this case there 
is one equilibrium - minimum E = Q that corresponds, for E = E™ = ^ln{gf^ g) 
to the value g{X, T) = g^ -constant. We have Y = c. 

2. Two-phase material (material that can exist in two stable phases). In this 
example function F{E) has two (locally) stable states g = go,gi, or E ^ 0, E = 

F{E) = Fo + ^CiE^iE + Q-^f + ^C2E^. (12.12) 

Then, for E = E"\ 

FiE) = Fo + lc,ilH9o'9)nlH9i'9)? + lc,{Unig^' g)^ . (12.13) 



We have 



F"{E) =C2 + ^Ci{{E + Q-^f + AE{E + Q-^) + E^). 



Young module in the state g — g^ is equal to 

1 _2 

Y(3^ C2 + -CiQ , 

while in the second stable state g — gi, 

n = C2 + ic4((Q + Q-^f + AQ{Q + Q-i) + Q^) = Yo + 3c4(l + Q") 

Thus, in a case of a wave of phase transition going along the bar, Young module 
changes by the amount Yi — Yq = 3c4(l + Q^). 



13. Appendix II. Variations 

Variations of some expressions for the Lagrangian (5.4-5.5) are calculated here 
and presented in a table. All terms in Lagrangian density L{G, 4>) will be refereed to 
the mass form dM — podcV = pQS^y\g\d'^X . In other terms we calculate variation 
/ f{A)dM by A. Resuh of variation has the form VdM : Sf{A)dM = VdM. 
In the calculations we repeatedly using the following standard relation (Jy^jyj = 



■^^^—g^'''Sgij (see, for instance, |H]). In the table below we present tensors V for 
different /. 

As an example of such a calculation we provide calculation of S [diVg{N)dM] 
where diVg{N) = -j^^{^\N'): 
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^./^. 



SdiVg{N)dM = diVg{N)po{^\g\SS+S^g''Sgij)+5i 



\9 



dXi 



{^\N'))dM 



d 



(diVg{N)S-'SS-—-jln{poS)SN' 



dX' 



d^Vg{N)l ~ \{diVg{N) ~ In' J^lnipoS)) 



{divg{N)S-^5S 



dXi 



ln{paS)5N' 



2 2 



2" dX' 

g^''SgAB)dM 

(13.1) 



g''"6gAB)dM 



since 



5{ 



■j{^\N'))dM = -{-^ 



-\{d^v,{N) - ^N'^lnipoS))g^^Sg 



9\J^{V\9\N'))dM-j^{poS)6{y^\N')d^X 



dXi 



dXi 



d 



AB 



dX 



jln{poS)5N']dM (13.2) 



Formula of variation of hR{g) y/\g\ in the 5th row of the table is taken from |llj . 
Prop.3.2. 

To find variation of the strain energy density f{G,E'^'^)poS^/\g\d'^X we first 
take variation of dM = poS^/\g\d'^X to get the first two terms in the last row of 
the Table, then - explicit variation by g if the strain energy function / depends 
on g not just through E^K Finally for variation by g through the strain tensor 

E'^'j = i5'^(C(</))Kj - OKj) we have 



SfiE^') 



df 
dE< 



SE^ 



1 df 



2qe<^i 



j5g'''C{ct>)Kj - 1 ^^ r .ia.kb 



"^dE^ij 

1 df 



j[-g'^g''^SgAB]C{^)Kj 



-lT^ii9'^9''''+9'''g''^)S9AB]Cicl^)KjSgAB = -jj^j^^^iKgiN + S'hgiM) x 
X {{g'W' + g'^g^^)SgAB)]C(c^)KjSgAB - -^( ^^ -^^ 



dE'^^ 



Of ^KA_ 



JB 



df KB , 9f ^i^A 



A^dE^'ijA 
g'''')Ci^)KjSgAB = -Is^^^'^SgAB, (13.3) 



where S**^"^^' is the symmetrization of the second Piola-Kirchoff Tensor 5* (see |6,), 
the last equality is proved in |22| . 



14. Conclusion 

In this work, we consider the intrinsic material metric tensor to be an additional 
parameter of state, i.e., an internal variable that characterizes material degradation 
and aging. The material metric tensor G is a conjugate (with respect to a particular 
Lagrangian) to the canonical Energy-Momentum Tensor (or to the Eshelby energy- 
stress tensor to some degree). 

Equations of metric evolution, (i.e., the aging equations), are derived as the 
Euler-Lagrange equation of a corresponding variational problem. Canonical energy- 
momentum tensor (or Eshelby Tensor) play a role of the source of metric evolution. 
This represents an alternative approach to numerous phenomenological damage 
models, which usually have more adjustable parameters than practical testing is 
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Term 


Variation by S, N^ , gjj 


/(l.9l) 


S-'f{\g\)5S + [ifM)g'' + f{\g\)\g\g'']Sgij 


fiS) 


[f{S) + S-'.f(S)]SS+if{S)9'-'Sgjj 


ll^ll^a 


\\N\\lS-'SS + 2NjSN' + [N^NJ + UN\\l9"]S9ij 


xiK'j) 




hR{g) 


S-'hR{g)5S + [-h{mc{g) - ^R{g)g) + -l^[^g{p^hS)g^" + Hess{pohS)]]5gAB 


dlVg{N) 


diVg{N)S-HS - ^ln{p^S)5N^ + \^lN^ ^ln{p^S))\ g'Hgjj 


fiG,E^') 


S-\fSS + [i/g^^ + ^e.p - ^S^^^^SsAB 



Table 2. Table of Variations. 



able to determine. Thus it is difBcult to validate the models since they can almost 
always be adjusted to reach an agreement with the experiment. In contrast, a 
variational approach prescribes a functional form of the aging equations, limits the 
number of constants (adjustable parameters) employed in the Lagrangian, provides 
a simple physical interpretation of the constants, and admits an essential experi- 
mental examination of the validity of the basic assumptions of the model. Particular 
examples (aging homogeneous rod, see Sec. 11 or [20], cold drawing (necking) [18], 
residual stress and others) can be analyzed theoretically and unambiguously tested 
in the experiments as a natural continuation of the present work. 

A natural development of this scheme requires the following: thermodynamical 
interpretation of the balance equation considered in section 10, especially the struc- 
tural entropy evolution manifested by the increase of the lapse function S{t) during 
aging; introduction of a temperature dependence of the material metric (based 
on an unpublished work by A.Chudnovsky and B.Kunin); development of models 
(" ground energy" F + kinetic potential A + possibly other metrical terms) charac- 
terizing a hierarchy of aging phenomena for specific materials; and especially, the 
development of models of phase transition (front propagation, fractal restructuring 
of materials, etc.). 

We would like to express our gratitude to Professor M. Francaviglia for his at- 
tention to this work. We would also like to thank Professor R.Tucker and the 
participants of his seminar in the Physics Department of Lancaster University, UK 
for the useful discussion. 
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